REGULAR SINGULAR STURM-LIOUVILLE OPERATORS AND 
THEIR ZETA-DETERMINANTS 

MATTHIAS LESCH AND BORIS VERTMAN 

Abstract. We consider Sturm-Liouville operators on the line segment [0, 1] with 
f^ ■ general regular singular potentials and separated boundary conditions. We establish 

Cn I existence and a formula for the associated zeta-determinant in terms of the Wronski- 

^ i determinant of a fundamental system of solutions adapted to the boundary conditions. 

This generalizes the earlier work of the first author, treating general regular singular 
potentials but only the Dirichlet boundary conditions at the singular end, and the 
>^ ^ recent results by Kirsten-Loya-Park for general separated boundary conditions but 

^\j . only special regular singular potentials. 
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1. Introduction and formulation of the result 

In this paper we will investigate the zeta-determinant of Sturm-Liouville operators 
of the form 

dP z/2 - 1 /4 1 
H = --— + ^ + -V{x), Rez/>0,xG (0,1), 

the regularity assumptions on V will be minimal. 

Such operators are sometimes also called Bessel operators. H is the prototype of a 
differential expression with one regular singularity and hence it appears naturally in the 
classical theory of ordinary differential equations with regular singularities [CoLe55]. 
The physical relevance of H stems from the fact that it arises when separation of 
variables is used for the radial Schrodinger operator in Euclidean space. 

Quite recently there has been a lot of interest in the inverse spectral theory of iJ, see 
e.g. [KSTIOb], [KSTIOa], [Car97] and the references therein. 

Our motivation for looking at the zeta-determinant of H comes from geometry: Spec- 
tral geometry on manifolds with singularities has been initiated by Cheeger in his sem- 
inal papers [Che79b], [Che83]. Manifolds with conical singularities are an important 
case study for this general programme. Separation of variables for the Laplacian on a 
cone leads to an infinite sum of Bessel type operators like H above. Recently, there has 
been a revived interest in extending the celebrated Cheeger- Miiller Theorem [Che79a], 
[Mul78] on the equality of the analytic torsion and the Reidemeister torsion to mani- 
folds with conic singularities [DaHuIO], [dMHS09], [HaSpIO], [Ver09], [VerIO]. 

The separation of variables mentioned above leads naturally to the problem of de- 
termining the zeta-determinant of a single regular singular Sturm-Liouville operator 
on the line segment [0, 1] with separated boundary conditions. We only make minimal 
regularity assumptions on the potential. Nevertheless, we establish existence and a for- 
mula for the associated zeta-determinant in terms of the Wronskian of a fundamental 
system of solutions adapted to the boundary conditions, see Theorem 1.5 below. 

We emphasize that for the calculation of the analytic torsion or the zeta-determinant 
on a cone additional considerations are necessary. This is because on a cone one has to 
deal with an infinite direct sum of operators like H. 

The fundamental results of Briining-Seeley in [BrSe85], [BrSe87] guarantee the ex- 
istence of zeta-determinants for regular singular Sturm-Liouville operators with Dirich- 
let boundary conditions at the singularity. However, loc. cit. require the potential to 
be of the form a{x)jx^ with a{x) smooth up to 0. For such operators Theorem 1.5 
was proved in [Les98] by the first author, generalizing earlier results by Burghelea, 
Friedlander and Kappeler [BFK95] to the regular singular setting. 

The method of [Les98] is limited to the Friedrichs extension at the singularity. In 
a recent series of papers Kirsten, Loya, and Park [KLP08b], [KLP08a], [KLP06] 
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were able to calculate the zeta-determinant for an explicit example of a regular singular 
Sturm-Liouville operator with general self-adjoint boundary conditions; cf. also the sub- 
sequent discussion by the second author in [Ver09] and in the appendix to [KLPOSa]. 
Their method, however, is based on an intricate analysis of Bessel functions and is 
therefore limited to their explicit potential. 

The main result of this paper combines and generalizes these two results, however 
only for scalar valued potentials. Since we deal with a rather general class of potentials 
it is natural that our method is closer to that of [Les98]. Special functions are used 
in this paper only implicitly as we are using the formula from [Les98] for the zeta- 
determinant of the Friedrichs extension of the regular singular model operator l^ = 

-& + (^' - 1/4) /o:^. 

The paper is organized as follows. In the remainder of this section we will in- 
troduce some notation, explain the basic concepts of regularized integrals and zeta- 
determinants, and we will formulate our main result. In Section 2 we derive the as- 
ymptotic behavior of a fundamental system for if, slightly generalizing a result due to 
Bocher [BocOO]. 

In Section 3 we study the maximal domain of H and its closed extensions with sep- 
arated boundary conditions. Let H{6q, Oi) be such an extension, 6q, 6i stand for the 
boundary conditions at 0, 1 resp. We give criteria under which it is possible to factor- 
ize H{9o,9i) into a product -D1-D2 of closed extensions of first order regular singular 
differential operators. We prove a comparison result for the Wronskians of normalized 
fundamental solutions for D1D2 and D2D1. 

In Section 4 we discuss the asymptotic expansion of the resolvent trace. We start 
with the Friedrichs extension. The resolvent of the Friedrichs extension Li, of the 
regular singular model operator l^ is explicitly known and rather well-behaved with 
respect to perturbations of the form X~^V. From [BrSe85] we only use the result 
that the resolvent of the model operator L^, has a complete asymptotic expansion. 
The expansion of the resolvent of L^, + X~^V then follows by a perturbation analysis. 
Boundary conditions other than the Friedrichs extension at are more subtle since the 
resolvent does not absorb high enough negative powers of x. For the resolvent of general 
boundary conditions we therefore use the factorization results of Section 3. In addition 
one needs to treat compactly supported L^-perturbations of factorizable operators. For 
this we employ a standard method of pasting together local resolvents, cf. [LMP09, 
Appendix] . 

In Section 5 we derive a variational formula for the dependence of the zeta- 
determinant under variation of the potential. The method is well-known [BFK95], 
[Les98]. However, due to the low regularity assumptions on the potential and due 
to the singularity of the operator the analysis becomes a little delicate. In particular 
we have to analyze the dependence of a normalized fundamental system (and its as- 
ymptotic behavior near 0) on the parameter. At the end of Section 5 we compile the 
established results to a proof of the main Theorem 1.5. 
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1.1. Function and Distribution Spaces. Following the requests of several of the 
referees we are going to specify in detail the notation for function and distribution 
spaces used throughout the paper. 

Let / C M be an interval, which may be of any of the possible forms (a, b), [a, b), (a, b] 
or [a, b] for real numbers a < b. Let 1° = I \{a,b} denote the interior of I. 

For a map f : I ^ E into some vector space E the support of /, denoted by supp /, 
is defined as the closure in / of {x G / | /(x) 7^ 0} 

-/ 



snppf:={xel\fix)y^0} . (1.1) 

supp / is always closed in / but not necessarily compact, since / might be non-compact 
itself. 

For spaces of continuous, respectively differentiable complex- valued functions we use 
the standard notation C(/), C'=(/), C°°(/), cf. e.g. [Hor90, Sec. 1.1, 1.2]. The space 
Cq{I),0 < k < 00, denotes the subspace of those / G C^{I) with compact support. 

The space C^(/°) carries a natural locally convex topology and its dual space &\I°) 
is called the space of distributions on 1° . 

For T G ^'(/°) one can define suppT [Hor90, Sec. 2.2]. For an arbitrary subset 
A C M one now writes ff' (A) = {T e ^'(M) | suppT C A], cf. [Hor90, Sec. 2.3]. 
For the half open interval / = (a, 6], e.g., T G S"{{a,h\) if there is a 5 > such that 
suppT C (a -|- 5, 6]. 

For distributions it also makes sense to talk about restrictions. If J C / are intervals 
and T G ^'(/°), we put T\j := T \ C^{J°). 

Let F be a map which assigns to each interval / C M a subspace E{I) C &'{I°). 
Furthermore, assume that F is compatible with restrictions in the following sense: if 
J d I are intervals and / G F{I), then f\j G F{J). Then Fioc(/) denotes the space 
of T G ^'{1°) such that T\k G F{K) for each compact interval K G L Furthermore, 

Fcomp(/):=Floc(/)n^'(/). 

Example 1.1. For an interval / C M we denote by L^(/), 1 < p < 00, the Banach space 
of p-summable (equivalence classes modulo equality almost everywhere) functions with 
respect to Lebesque measure; for / G L^il) the norm is given by 

\LP ■= \ I \f{x)\Pdx^ 



LP{I) is naturally embedded into S>'{I°) by identifying / G F^^I) with the distribution 

C-(r)9 0^ /"Z-^; (1.2) 



J I 
needless to say that (1.2) is independent of the choice of a function representative of the 
class /. The support of / G U^I) is now defined as the closure in I of the support of 
the corresponding distribution in ^'{1°). For continuous functions C{I) C L^il) (each 
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L^-class has at most one continuous representative) the latter definition of support 
coincides with (1.1), assuming that / is closed. 

The assignment I i— )■ U'{I) is an example for the map F discussed above. Hence, 
-^comp(-^) ^^d -Z^foc(-^) ^^^ defined. Note that although L^^I) = L^^I), we only have 

^om^in C ^.^^(I), LU7) C LUI). 

Sobolev spaces will only be used in the Hilbert space setting p = 2. We write 
H^{I) for the Sobolev space W^^^{I) of those / G L'^{I) C ^'(/°), for which all weak 
distributional derivatives d^ f, 1 < j < k, taken a priori in ^'(/°), are actually in L'^{I). 

1.1.1. The Schatten Ideals. For a Hilbert space "H we denote by i3('H) the space of 
bounded and by /C('H) the space of compact operators on "H. For 1 < p < oo let 
B^iTi) C /C('H) be the von Neumann-Schatten ideal of p-summable operators, cf. e.g. 
[Ped89, Sec. 3.4]. For T G B'^{'H) the p-norm is given by 

\\T\\, := (Tr(T*T)^/2) '^^ = ( ^ \^'^) "'. 

AGspcc T*T 

Tr denotes the trace [Ped89, Sec. 3.4]. We will only need p = 1 and p = 2. Operators in 
B^{l-i) are called trace class operators and elements of B^{'H) are called Hibert-Schmidt 
operators. To avoid possible confusion with the L^-norm of functions, we write || ■ ||tr 
for the trace norm || ■ ||i and || ■ ||hs for the Hilbert-Schmidt norm || • ||2. 

1.1.2. Regularized Integrals. Let us briefiy recall the partie finie regularization, cf. 
[Les97, Sec. 2.1], [Les98] and [LeTo98], of integrals on M+ := [0,oo). Let 
/ : (0, oo) — !■ C be a locally integrable function. Assume furthermore, that for x > Xi 
we have a representation 

N 

f{x) = j2fr^"'+9{^)^ (1-3) 

with real numbers aj, numbered in descending order with a^v = —1, and g G L^[xi, oo). 
Then 

R 

f{x)dx 

XI 

N-l .oo .oo (1-4) 

=:^—^R'''+^ + f^logR+4 f{x)dx + o{l), asR^oo. 

o{l),R — )■ cxD, is the usual Landau notation for a function of R whose limit as i? — > oo 
is zero; here we have explicitly o(l) = /„ g. The regularized integral ^ f{x)dx is 

therefore defined as the constant term in the asymptotic expansion of / f{x)dx as 
i? ^ oo. 

If for < X < Xq we have a representation 

M 

f{x) = Y.fy^ + h{x), (1.5) 



MATTHIAS LESCH AND BORIS VERTMAN 















• 


• 








• 


• • 








y . . 


.\ 






• 




• 


r 




• 









Figure 1. The contour of integration F. 



with real numbers /3i < (32 < ■ ■ ■ < Pm = —1, and h e -^-"^[O, xq], then 



Xo 



f{x)dx 



M-l rO 



y2ir^^^'^'-fM^og6 + 4 f{x)dx + o{l), as (5^0, 

7 = 1 Pj + -*- Jo 



;i.6) 



and the regularized integral j^ " f{x)dx is defined as the constant term in the asymptotic 
expansion of J^ f{x)dxa.s6-^0. 

Now assume that / satisfies (1.3) and (1.5). Since / is locally integrable, it is clear 
that (1.4) holds for any Xi > and (1.6) holds for any xq > 0. One then puts for any 
c> 

/■OO /-C /"OO 

(1.7) 



f{x)dx:=-j- f{x)dx+-j- f{x)dx, 
Jo Jc 

and in fact the right hand side is independent of c > 0. 

1.2. The zeta-determinant. Let i7 be a closed not necessarily self-adjoint operator 
acting on some Hilbert space with spec(— if) fl R+ finite, ^ specif. We assume that 
the resolvent of H is trace class, and that for z E M., z > zq > max(spec(— if ) fl M+) 



Tr(ii + z 



-1 



'2 ^ T-./ \ 

+ - + i? 02 

z z 



with 



lim zR{z) = 0, 

\R{z)\dz < OO. 



2— >00 
fOO 



(1.8) 

(1.9) 
1.10) 



^0 



Let F be the contour as sketched in Figure 1, where the bullets indicate the eigenvalues 
of —H. Fix a branch of the logarithm in the simply connected domain C \ {— F}. Note 
that the previous definition (1.7) of the regularized integral can easily be adapted to 
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functions defined on the contour F, since there are < xq < xi < cxd such that [0, xq] 
and [xi, oo) are contained in (the image of) F. Consider for fixed s G C the function 

f^{x) ■=x" Tr(/7 + x)-\ 

In view of (1.8) and (1.10) it satisfies (1.3) if Re s > 0. Furthermore, since H is assumed 
to be invertible, the function x H- Ti{H + x)~^ is smooth up to a: = and its Taylor 
expansion at x = shows that /^ satisfies (1.5) for all s G C 

Exploiting the definition of jf. it is now not hard to see, cf. [LeTo98, (2.30)], that 
for 1 < Re s < 2 the zeta-function is given by 

Ch(s):= y \-' =^-^^^^4 x-'TiiH + xy^dx. (1.11) 

Aespcc//\{0} ^ 

Furthermore using the asymptotic expansions as x — )■ oo and x — > of x~* Tr(i7 + 
x)~^ one deduces that the right hand side of (1.11) extends meromorphically to the 
half plane Res > 0, [LeTo98, Prop. 2.1.2]. The identity (1.11) persists except for the 
poles of the function s t-)- ^j^C-ffls). Thanks to (1.10) the function C,h is differentiable 
from the right at s = and one puts 

logdetf if := -Ch(0) = - -[ Tt{H + x)-^dx. (1.12) 

det(^ H is called the zeta-regularized determinant of H. For non-invertible H one puts 
deti^H = 0. With this setting the function z i— )■ det(^{H + z) is an entire holomorphic 
function with zeroes exactly at the eigenvalues of —H. The multiplicity of a zero z 
equals the algebraic multiplicity of the eigenvalue z. 

1.3. A regular singular operator. We now introduce the class of operators we are 
going to study in this paper. Put 

rf2 z/2 - 1 /4 
/. = -— + ^^, Re^>0, (1.13) 

acting as an operator in the Hilbert space L'^[0, 1], a priori with domain C^{0, 1). We 
will study perturbations of l,y of the form 

H = l, + X-^V, (1.14) 

with suitable conditions on the operator V to be specified below. X denotes the function 
X{x) = X. We view V resp. X~^V as a multiplication operator on functions on the 
unit interval. In order not to overburden the notation we will in general not distinguish 
between a function / and the operator of multiplication by /. 

Definition 1.2. 1. For an interval / C M we denote by AC^{I),k > 1, the space 
of / G C{I) C ^'(/°) such that d^ f G C(/),0 < j < k - I, and d'^f G Aoc(^)- 
AC^{r) = AC{I) is the well-known space of absolutely continuous functions. Note 
that for this definition it matters whether a boundary point p G (9/ belongs to / or not. 
2. Denote by % the space of those V G Ly^ci^, 1) such that 

V-\og{j)eL%l], ifz/^0, (1.15) 

r-log'(2) eL'[0,l], if 1^ = 0. (1.16) 
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A natural norm on "^ is given by 










\V\r.= log( 


2)y .., 


if 1/ ^ 0, 




(1.17) 


\V\-K= log( 


2> ^ Li' 


if z/ = 0. 




(1.18) 


■^ is a Frechet space with seminorms | • 


|r. and 


I^[l/n,l-l/n]l^^'^ = 


= 2,3,.., 





3. Finally, let ^ be the space of those / G AC^{Q,l) such that f\Xf" G 

Some of the results will hold under the weaker hypothesis V G log(2)X' ^ X-- Unless 
said otherwise, function spaces consist of complex valued functions. 

In Section 2 we will prove the following refinement of a Theorem of Bocher [BocOO] 
(Theorem 2.1 and Proposition 2.6): 

Theorem 1.3. Let V G Yjieu.H = l^ + X'^V , Rez/ > 0, and let ui = u + 1/2, z/2 = 
— Z/ + 1/2 be the characteristic exponents of the regular singular point of the differential 
equation Hg = 0. Then there is a fundamental system gi,g2 of solutions to the equation 
Hg = such that 

g,{x) =x''^giix), (1.19) 

,,(x)=|-^^"^^f^)' '^'^°' (1.20) 

\^^/x\og{x)g2{x), if z/ = 0, 

where gj ^ s^ . 

The spectra and fundamental system of solutions to Bessel type Sturm-Liouville dif- 
ferential expressions on finite intervals have also been studied (mainly in connection with 
the inverse spectral problem) in a number of recent publications [AHM07], [Car97], 
[EvKa07], [KSTIOa], [KSTIOb] and [Ser07]. 

1.4. Separated boundary conditions. Denote by Hq the differential expression H 
restricted to C^(0, 1) C -^^[0,1]. Let H^ be the formal adjoint of Hq. This is the 

differential expression — ^ + ^ ~2 — h ^^i^) acting on C^(0, 1). Hq is symmetric if 
both z/ G M and V is real valued. 

As usual we denote by -ffmin = Hq the closure of Hq and by H^aax = (Hq)* = 
(if^in) . For convenience we introduce the left minimal and right maximal domain 
!^l{H) as the domain of the closure oi H \ C^{0, 1]. The left maximal and right minimal 
domain ^r{H) is defined accordingly with C^{0,1] replaced by {/ G i^(ifniax) | 
supp/ C [0, 1) compact}. Note that by the definition of support in (1.1), compactness 
of supp / C [0, 1) means that supp / has a positive distance from the point x = 1. 

Although there is no simple Weyl alternative in the non-self-adjoint context, we say 
that X = (resp. x = 1) is in the limit point case for H if ^l{H) = ^{H^g,^) (resp. 
&r{H) = ,^(iJmax))- Otherwise, we say that it is in the limit circle case. 

We will see in Section 3 that there are continuous linear functionals Cj,j = 1,2, on 

^(il^max) such that for / G ^(i/max) 

/ = Ci(/)(7i + C2(/)^2 + 7, fe^LiH), (1.21) 
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where gi,g2 are defined in (1.19) and (1.20). 

A boundary condition at the left end point is therefore of the form 

Bo^ef ■■= sine- ci{f) + cose- C2{f) = 0, < 9 < tt. (1.22) 

6 = gives the Dirichlet boundary condition (Friedrichs extension near 0). 

It should be noted here that is in the limit point case for H if and only if Re z/ > 1. 
In this case g2 ^ L^[0, 1], C2 = 0, and hence &l{H) = ^(ifmax). Thus if Rez/ > 1 we 
consider only the case 6 = 0. Boundary conditions such as (1.22) at the singular end 
point have been studied in depth by Rellich [Rel44] and extended by Bulla- Gesztesy 
[BuGe85]. 

From the well-known fact that a linear second order ODE with L^-coefficients has 
y4C^-solutions it follows in view of our assumptions on V that ^{H^i^^) C AC^{0, 1] 
and hence 1 is always in the limit circle case for H. At the right end-point we therefore 
impose boundary conditions of the form 

Si,0/:=sin^-/'(l) + cos^-/(l) = 0, < ^ < vr. (1.23) 

For each admissible pair (^o, ^i) ^ [0, vr)^ (0 < 6'o < vr if Re z/ < 1, 6^0 = if Re z/ > 1) we 
obtain a closed realization H{6q, 6i) of the operator with separated boundary conditions 
Bog^jBig-^. All eigenvalues of H{6o,6i) are therefore simple. 

Under the technical assumption that V is of determinant class, see Definition 4.4, 
which is satisfied for all real valued potentials V^ G "^ we can prove (Theorem 4.3 and 
Theorem 4.10) 

Theorem 1.4. Let u > 0, V E 'fy, and assume that 6q = or that V is of determinant 
class and u > 0. Then the resolvent of H{6q,6i) is trace class. Moreover, there is 
a zq > such that H{6q, 6i) + z is invertible for z > zq and there is an asymptotic 
expansion 

-1 a h 

\/z z 
with 



Ti{H{6o,6i) + z) ^ = -= + - + R{z), z>ZQ,ze^, (1.24) 



lim zR{z) = 0, (1.25) 

\R{z)\dz <oo. (1.26) 



z— >oo 

r"00 



In view of this Theorem we may define det^ (if (6'o, ^i)) according to (1.12). 

1.5. The main result. To explain our result we need to introduce the notion of a 
normalized solution at one of the end points. First, we define an invariant of the 
boundary operator Bjq: 

fz/, if = 0, , , 

— z/, if < t^ < tt; 



resp. 



,1/2, if = 0, , , 

/xi:=/x 5i0 ■= { , 1-28 

^ ^^ '^ ^ -1/2, ifO<0<7r. ^ ^ 
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To explain the ±1/2 we note that the right end point may artificially be viewed as a 
regular singular point with v = 1/2. Hence fij depend in fact on 6 and the characteristic 
exponent of the regular singular point. 

A solution of the homogeneous differential equation Hg = is called normalized at 
with respect to the boundary operator 5o.e if Bo^gg = and if g{x) ~ x^°~^^^'^, as 
X — )■ 0; here we use the notation 

fix) 
f{x) ~ h{x), as x -^ xo :^ lim f^ = 1. (1.29) 

x^xo h[x) 

Similarly, g is called normalized at 1 with respect to the boundary operator Biq if 
Bi^g = and if g{x) ~ (1 — 2;)^i+^/^ as x — )■ 1. It is straightforward to check that there 
is always a unique normalized solution. 

Theorem 1.5. Let Bj£.,j = 0,1 be admissible boundary operators for H. Under the 
same assumptions as in Theorem 1.4 the zeta-regularized determinant of H{9q^9i) is 
given by 

detc(H (60,61)) = — — —- ^^— -Wi-^^ip). (1.30) 

fV V u, ijj 2«'+^'ir(//o + l)r(/ii + 1) ^^'^' ^ ' 

Here, if, tp are solutions to the homogeneous differential equation Hg = such that 
if is normalized for -Bo,6»o T^^ ^) '^'^^ V" ^■^ normalized for Biq^ (at 1). Furthermore, 

W{ip, (f) = ip(p' — il)'(p denotes the Wronskian of ip, (f. 

Theorem 1.5 has a relatively straightforward extension to the case where the potential 
has regular singularities at both end points. The proof does not require any essentially 
new idea; the details, however, are a bit tedious and are therefore left to the reader, cf. 
Remark 5.6. 

The case i/ = and 6q > 0, which is not covered by Theorem 1.5, requires specific 
analysis of unusual singular phenomena in the trace expansion of H, as observed first 
by Falomir, Muschietti, Pisani, and Seeley [FMPS03]; see also the nice elaboration by 
Kirsten, Loya and Park in [KLP06]. The discussion of the zeta-determinant in this 
case therefore requires another publication. 

To outline the proof of Theorem 1.5 we first observe that if Di, D2 are closed operators 
in a Hilbert space then spec D1D2 U {0} = spec 1^2-01 U {0} and, even more, nonzero 
eigenvalues of D1D2 and D2D1 have the same multiplicity. Hence if both D1D2 and 
D2D1 satisfy the general assumptions of Section 1.2 then for z E'C 

deii^{DiD2 + z) = z'^dei(^{D2Di + z), rf := dimker ADz - dimker Da^i- (1-31) 

We will show in Proposition 3.5 that H{6q, 61) can always be written in the form 

H{6o,6i) = DiD2 + W (1.32) 

with a compactly supported L^-potential W and Di,D2 suitable closed extensions of 
the operators 

d oj' , d u' ,_, „„, 

rfi = — + -, ci2 = -— + - 1.33 

ax 00 ax 00 

with a certain function u which is singular at 0; its properties will be described in detail 
in the text. The crucial point is that for the interesting case ^o > one can choose 
Di,D2 in such a way that D2D1 also is an operator to which Theorem 1.4 applies 
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and such that the boundary condition at is the Friedrichs extension. The Friedrichs 
extension at is much better behaved and can be treated for our class of operators 
basically as in [Les98]. The proof is completed then by employing a variation result 
for the behavior of the zeta-determinant under variation of the potential W (Theorem 
5.4). 



2. The fundamental system of a regular singular equation - Bocher's 

Theorem 

Consider the following regular singular model operator 

£ z/2-1/4 

acting on C^(0, 1) C L'^[0, 1]. z/ is a complex number for which without loss of generality 
we may assume Rez/ > 0. 

We are interested in perturbations of the form 

H:=l, + X-W, (2.2) 

with V G -/^ioc(0' 1) ^^'1 -^ denoting the function X{x) = x. In this section we are 
concerned with the description of the asymptotic behavior as a; —> of a fundamental 
system of solutions to the equation Hf = 0. 

If V is analytic, then the classical theory of ordinary differential equations with regular 
singularities, cf. e.g. [CoLe55, Chap. 5], apphes and the characteristic exponents of 
the regular singular point at x = are given by 

z/i = z/ + l/2, z/2 = -i/+l/2. 

Furthermore, there is a fundamental system of solutions to Hf = of the form 

fi{x)=x'fi{x), f2{x) = < ^ / X 7/ X -r n ^^-^^ 

[Vxlog(x) /2(X), lfl/ = 0, 

where /j, j = 1, 2, are analytic functions with /j(0) = 1. The normalization of solutions 
is chosen so that 

W{hj2) = hn-f[f2 = l. (2.4) 

It is less known that already M. Bocher [BocOO] investigated regular singular points 
of ordinary differential equations with non-analytic coefficients. For Bessel operators 
with L^ potentials a thorough analysis of the fundamental system of solutions was made 
e.g. by Carlson [Car97]. Bocher's result reads as follows. 



+ -V{x) 

X 



Theorem 2.1. 


M. Bocher] Let 












H = - 


dx^ 


+ 


Z/2 


-1/4 
a;2 
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where z/ G C, Rei/ > 0, and V G \og[^^Yy. Then the differential equation Hg = has 
a fundamental system of solutions gi,g2, such that 




^x'^^M^), ifz/^0, 
xlog{x)g2{x), if z/ = 0, 



(2.5) 
(2.6) 



where g, G C[0, 1], g,{0) = 1 forj = 1, 2. 
Furthermore, 

g[ix) = uix^'-'h^ix), (2.7) 

, i-^^x^^-'h^ix), if.^0, 

l27^M^)^2(a:), ifz^ = 0, ^'-'^ 

w/iere /i^- G C[0, 1], hj{0) = 1 /or j = 1,2. 

Finally, with these normalizations 

W{gi,g2) = gig'2- 9'i92 = 1- 

Remark 2.2. In the case z/ = the theorem as stated is shghtly more general than 
[BocOO], where V^ log G L^[0, 1] is assumed. Moreover, note that the conditions on the 
potential V in the theorem are satisfied whenever y G L^[0, l],p > 1, or more generally 

Ve%. 

We briefiy sketch a proof of Theorem 2.1 in modern language. Being self-contained 
is not the only reason for presenting the proof in some detail: the method of proof will 
allow a more precise analysis of the regularity properties of gj (see Prop. 2.5 and Prop. 
2.6 below) which will be needed later on. Furthermore, the method will be needed for 
deriving the variation formula for the zeta-determinant under variation of the potential 
(Section 5). 

Proof. The regular singular operator ly has the following fundamental system of solu- 
tions to luf = 0: 

Lx''-^ if 7/ =^ 

X , 1 z.^U, 

Va;logx, it z/ = U. 

For the Wronskian we have W{fi, f2) = /1/2 — /1/2 = 1- For a solution to Hg = 
we make the ansatz 

g^{x) = fi{x) + x'''(j){x) = x''i(l + 0(x)). 

Plugging this ansatz into the ordinary differential equation Hg = yields for iIj{x) = 

-r{x) + ^^^^^i;{x) = -^V{x)[f,{x)+ij{x)], (2.9) 
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thus 

0(a;) =-x-^i/i(a;) [^ Uy)-Viy)y-'^[l + <P{y)]dy 

Jo y 

+ x-^72(x) r h{y)-V{y)y''^[l + <P{y)]dy ^^'^^^ 

Jo y 

=:{K,Vl){x) + {K,V<l)){x), 
where K,y is the Volterra operator with the kernel 

fc,(x,i/)=i^(l-x-V^), y<x,ifu^O, (2.11) 

koix,y) = -\og{y) + \og{x), y<x,iii^ = 0. (2.12) 

We view K^V as an operator on the Banach space C[0, 1]. Indeed, for any (p G C[0, 1] 
one easily checks 

\KV<Pix)\<^ r \V{y)My)\dy, 
W\ Jo 

\KoV<Pix)\ < 2 r \logiy)\\Viy)My)\dy, 
Jo 

From (2.13) and (2.14) one infers by induction 

\iKVr<t>ix)\ < Y^mU[o J r \V{y)\dy 
|z/|"n! \Jq 

on , px 

\{KoVr<P{x)\ < -||0||^,[o,x](y \V{y)logy\dy'^ 

Hence for any u G C, Re z/ > 0, the Volterra operator K^^V is a bounded operator 
on C[0, 1] with spectral radius zero. Consequently the equation (2.10) has a unique 
solution in C[0, 1] given by 

= (/ - K^Vy^K^Vl. (2.17) 

Moreover, by (2.15) and (2.16) one has 

[C2J0 \V{y)\ogy\dy, if z/ = 0, 

for some constants Ci,C2, not depending on V. This proves that 

gi{x) = x"^ {1 + (f>{x)) = x^^gi^x), 

is indeed a non-trivial solution to Hg = with gi G C[0, 1] and 5^i(0) = 1. To see (2.7), 
note that by (2.10) (p is absolutely continuous in (0, 1) with its derivative given by 

ex 

'2u-l I „,2u^ 



if z/ ^ 0, 


(2.13) 


if z/ = 0. 


(2.14) 


if z/ ^ 0, 


(2.15) 


if z/ = 0. 


(2.16) 



<f)\x) = x-'""-' / y"'V{y){l + (p{y))dy, for all Rez/ > 0. (2.19) 

Jo 

This implies 

W{x)\<- r\V{y)\dy 
X Jo 

and thus (2.7) and the claims about gi are proved. 
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The second solution g2 can now be constructed as usual by putting near x = 

g2{x) = C{x)g,{x), (2.20) 

where 



X 



\og{x)C{x), if z/ = 0. 



C is continuous over [0,Xo) and Xq G (0, 1] is chosen so that gi{y) 7^ for < y < Xq. 
Such an xq exists, since 5^i(0) = 1 and gi G C[0, 1]. It is then straightforward to check 
that g2 extends to a solution to Ig = on (0, 1] which has the claimed properties. D 

Now we come to the aforementioned improvement of the regularity properties oigj{x) 
as X — ;• 0. 

Lemma 2.3. Let f G AC^{0,1) with f',Xf" G L^[0,l]- Then (X/')(0) = 0. This 
holds in particular for / G ^2/ (of. Def. 1.2j. 

Proof. By assumption the function F{x) := xf'{x) — J^ h{s)ds,h := (X/')', is locally 
absolutely continuous and F' = 0, hence F = (X/')(0) =: c is constant. Thus 

f(x) = - + - r h{s)ds. (2.22) 

X 3; Jo 

By assumption we have /', h G -^^[0, 1]. Thus 
/(I) - fix) = [ f'{s)ds 



r [^ (2.23) 

= — clog X — log X / h{s)ds — i h{s) log sds 
Jo Jx 

= — clogx + o(logx), as X — )• 0, 

since for < 5 < 1 we have | J^ h{s) log sds\ < Cs + \ logx| Jg \h\. Because the left hand 
side of (2.23) is bounded it follows that c = 0. 

The last claim follows, since / G ^ implies /', {XfJ G L^[0, 1]. D 

Lemma 2.4. Let a G C and p G log(2)'^ea- Put 

{X 
(2.24) 
-x-""-^ Jy''p{y)dy, if Re(a) < 0. 
X 

Then we have 

feL%i]nACio,i), xf'eL'[o,i], (x/)(o) = o. 

//p e y^ea then f,Xf' G log{^y'L^[0, 1], that is JJ e ^. 
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Proof. Integration by parts shows easily that / G L^[0,1] (resp. / e log(2) -^^[0,1] 
if p e yjiea)- Moreover, clearly / and hence also Xf are both locally absolutely 
continuous in the interval (0, 1]. Furthermore, we have 



Xf = (Xf)' -/ = -(« + 1)/ + p e L% 1] 

(resp. G \og{jy'L'[0, 1] if p G fkea). 

(X/)(0) = follows from Lemma 2.3 applied to J^ f. D 

Proposition 2.5. In the setup and notation of Theorem 2.1 we have for j = 1,2 

fi?, G AC[0, 1], X^. G AC[0, 1], (X^.)(0) = 0. 

Proof. We have for the first fundamental solution 

gi{x) = x^^gi = x^^i^l + 0(x)), 

where by is given by (2.10). The claim about gi now follows from Lemma 2.4 and 
the explicit form of the derivative (2.19). To prove the claim for g2, recall that for some 
Xq G (0,1], such that gi{y) ^ ioi < y < Xq, the second fundamental solution is 
given by 

92ix) = -^i(x) / giiyy^dy. 

J X 

If 1/ ^ 0, then 

g2ix) = -2vx''~^''^g2{x) 

r» (2.25) 

= 2ug,{x)x"' y-"'-%{y)-'dy=:U^)f{x). 

J X 

In view of the statement being proved for gi before and since the claimed properties 
are preserved under multiplication, it suffices to prove the claim for /. Integration by 
parts gives 

fix) = -x'^ iy'"'giiyr')C +X'-' p y''^piy)dy 

i>xo 

= c{xo)x'^ + g,{x)-' + x"" / y-"'piy)dy, (2.26) 

J X 

where p = (gi'^)' G L-^[0,Xo]. 

The first two summands are a priori in ylC[0, xq]. The last one is /1C[0, xq] by Lemma 
2.4. Furthermore, from the definition we infer 

Xf' = 2uf-2ugi'eAC[0,xo]. 

(X/')(0) = then follows from Lemma 2.3. 

Clearly, g2 and Xg'2 are locally absolutely continuous in the whole interval (0, 1] and 
hence the claim is proved for g2 and u ^ 0. 
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Finally, for i/ = we have 

y/x log X J^ 

-gii^) y~%iyy^dy=:gi{x)f{x). 



logx^ 
Again, it suffices to prove the claim for /. We compute with p := (gi'^)' G L^[0,xo] 



as before 



1 r^o 

f{x) = - 1 — / y'%iyy^dy 

log X J^ 

= - 1^ l°g(^)5^i(^)"ir + i^ I ^ogiy)piy)dy 
-"^"^"^ + Mx)-' + T^ r ^og{y)p{y)dy. 



log X log X 

From this one checks that / G AC[0,Xo] and hence g2 G AC[0, 1]. Furthermore 

xf{x) = -f^ - ■ 2/ X / ^og{y)p{y)dy, 
log X log (x) Jx 

and differentiating this again shows {Xf')' G i^^[0, 1]. The remaining claims now follow 
as in the case z/ 7^ 0. D 

Finally we prove the following refinement of the properties oigj,j = 1,2, which will 
be crucial for the rest of the paper. 

Proposition 2.6. Let V G i^Rcu- Then, in the notation of Theorem 2.1 we have for 
J = 1,2 

^.log(^),X^'log(2)GLM0,l], 

i.e. gj & £/ . 

Proof. We prove the result only for u ^ and leave the case i/ = to the reader. The 
result for z/ = will not be used in the rest of the paper. 

Recall that gi{x) = x'^^gi{x) = x^^^l + (/)(x)), where G y4C[0, 1] (Prop. 2.5) is given 
by (2.10) and observe that by (2.19) we have 

-2i/-l I „,2v 



<P'{x)=x-"'-'J y"'r{y)dy, (2.27) 



with r := V ■ (1 + 0), rlog(2) G -^^[0, 1] since V G '^ci/- The claims about gi now follow 
from (2.27) and Lemma 2.4. 

Recall from (2.25) g2{x) = gi{x)f{x) with / given by (2.26). Differentiating the 
latter we find 

l>XO 

fix) = c(xo)x2-i + 2Z/X2-1 / y-"^ig-'yiy)dy. (2.28) 



The first summand is still in L^[0,Xo] after multiplying by log (2)- For the sec- 
ond summand note that (gi^yiog[-) = — 25^^^5f^^ log(2) G L^[0,a;o] by the inclusion 
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g[\og[-) G L^[0, 1], proved above. Hence we can apply Lemma 2.4 to the second 
summand to conclude /'log(2) € L^[0,a;o]. Differentiating (2.28) we infer similarly 
X/"log(2) G L^[0,a;o]. Then one easily checks the claimed properties for the product 
92 = gif- Since g2log(^-^) and Xgl^logi^-^) are locally integrable in the interval (0, 1] we 
reach the conclusion. D 



3. The maximal domain of regular singular operators 

We continue in the notation of the preceding section and consider the regular sin- 
gular Sturm-Liouville operator H with the fundamental system {gi,g2) of solutions to 
the differential equation Hg = (cf. Theorem 2.1). We will freely use the notation 
introduced in Section 1.4. 

We have the following characterization of the maximal domain of H, compare 
[BrSe88] and [Che79b] and the basic discussion of the second author in [Ver09, 
Proposition 2.10]. Note that it holds under a slightly weaker assumption on the poten- 
tial (y G log (2)"^) than the one imposed in the rest of the paper. 

Theorem 3.1. Let ly he the operator (2.1) and let H = l^ + X'^V with V G \og[^)Yy, 
Rez/ > 0, and let gi,g2 be the fundamental system to Hg = of Theorem 2.1. Let f be 
a solution of the ordinary differential equation 

Hf = -f" + qf = geL'[0,l], (3.1) 

q{x) = '^^^^ + -Vix). 

Then f e AC\0,1] and 

fix) = c,{f)g,{x) + C2{f)g2{x) + /(x), (3.2) 

for some constants Cj{f),j = 1,2, depending only of f , 

fix) = O (a;^/2 log(a;)) , f'{x) = O [x^''^ log(x)) , a; ^ + . (3.3) 

Remark 3.2. We emphasize that the solution g\ is completely determined by the 
equation (2.5) and therefore canonical. However this is not so for gi- Surely, any 
function g2 + \g\ also satisfies (2.6). We mention this because as a consequence the 
functional C2 (!) is canonically given while C\ depends on the choice of g2- 

Proof. We first note that it is well-known that solutions to linear differential equations 
with Ll^^ coefficients are locally absolutely continuous. Therefore a solution / to (3.1) 
is absolutely continuous in the interval (0, 1] and from f" = g~ qf one then infers that 
/' is also absolutely continuous in (0, 1]. 
For xo G {0, 1} 

f{x)=gi{x) g2iy)g{y)dy - g2ix) gi{y)g{y)dy (3.4) 

Jxo Jo 
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is a solution to (3.1); note W{gi,g2) = 1- Depending on z/ we will choose xq such that 
(3.3) is satisfied. We first note that by applying the Cauchy-Schwarz inequality 

\92ix) gi{y)g{y)dy\<\g2{x)\n \gi{y)\'^dyj \\g\\L2 

'0(x3/2), ifz/^0, 

0{x^lHog{x)), ifz/ = 0, ^ 



(3.5) 



Furthermore, if Re z/ > 1, we put xq = 1 and find 

\9i{x 



C 92{y)9{y)dy < Cx^'^^+'l^lf y~^'''^'^^Hy)"^\\gU. 

Jx Jx 

0(a;3/2|log(a;)|i/2), ifRez/ = l, 



(3.6) 



0(a;3/2), ifRez/>l, 



Finally, if < Re z/ < 1, we put xq = and estimate 



92{y)9{y)dy\ < Cx'^'^^^^l^^j^ \g,{y)\'dyf"\\g\\L. 

0(x=^/Mog(x)), ifz^ = 0. 
This proves the estimates for /(x). Differentiating (3.4) we find 



(3.7) 



f\x)=9'i{x) 92iy)giy)dy-g'2{x) gi{y)g{y)dy, 

Jxo Jo 

and (2.7), (2.8) together with (3. 5), (3. 6) and (3.7) immediately give the claimed esti- 
mate for /'. Thus / is a solution to (3.1) satisfying (3.3). (3.2) is now clear. D 

Remark 3.3. The above proof shows that for z/ 7^ 0,Rez/ 7^ 1, the estimate (3.3) can 
actually be replaced by 

fix) = 0{x^^^), fix) = 0(x^/'), X ^ 0+ , (3.8) 

and for Re z/ = 1 by 

fix) = 0(x=^/2| log(x)|i/'), fix) = 0(xi/'| log(x)|^/2)^ X ^ + . (3.9) 

Corollary 3.4. Under the assumptions of Theorem 3.1 there are continuous linear 
functional Cj,j = 1,2 on ^{H^i^^) such that for f G ^(ifmax) 

f = c,if)g^ + C2{f)g2 + l (3.10) 

with f G ^l{H) (cf. Section lA). Let c^j be the corresponding functional for Hq. Then 
we have for f G ^{H^^,),g G &{H, 



max 7 



(H^^J^g) - {f,Hl,^g) = -/'(l)^(l) + fil)g'il) + C2if)c\ig) - cr{f)c\{g). (3.11) 
Finally, is in the limit point case for H if and only z/Rez/ > 1. In this case C2 = 0. 
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Proof. (3.11) follows easily from (3.2), (3.3) and the Lagrange formula. The formu- 
las (3.11) and (3.3) show that / G ^l{H)', the latter was defined in Section 1.4. 
Now, it follows from (3.2) that the quotient space ^{H^^x) / ^l{H) is spanned by 
g^ + !^L{H),g2 + &l{H) if 92 G L^[0, 1] and by g, + ^l{H) if g^ ^ ^^[0, 1]. This im- 
plies the continuity of the functionals Ci, C2 on ^(i^max)- Finally, C2 = if and only if 
g2 ^ -Z^^[0, 1]. The latter is equivalent to Re z/ > 1 and the claim is proved. D 

As already outlined in Section 1.4 we obtain a closed extension of Hq with separated 
boundary conditions by choosing boundary operators Bjg. 

5o,0j:=sin^o-Ci(/) + cos^o-C2(/), 
BifiJ := smOi • / (1) + cosOi ■ /(I), 

Bqq^ depends on the choice of a fundamental system, cf. Remark 3.2. To treat the 
limit point and limit circle cases at in a unified way we call a pair of boundary 
operators admissible if 61 G [0, vr) and either (6*0 = and Re i/ > 1) or (6*0 G [0, vr) and 
< Rei/< 1). 

Given an admissible pair -Bj,ej, j = 0, 1, of boundary operators we denote by H{Oq, 61) 
the closed extension of Hq with domain 

i7(^o,^i):={/G^(i^max) I 5„e,/ = 0,j=0,l}. (3.13) 

If z/ G M and V is real valued then Hq is symmetric and H[6q,6i) is self-adjoint. If 
Re 1/ > 1 then all self-adjoint extensions are obtained in this way. If Re z/ < 1 there also 
exist self-adjoint extensions with non-separated boundary conditions. These extensions 
will not be studied in this paper. 

3.1. Factorizable operators. Next we investigate when H can be factorized as did2 
with dj = ±^ + ^. For simplicity we confine ourselves to the case u y^ 0. Clearly, with 
some modifications one has similar results for u = 0. 

We have seen in Proposition 2.6 that if V" G yReu and w is a solution to the differential 
equation Huj = then uj{x) = x^^^^'^oj{x) with to G £/ and fi = ±z/. 

Conversely, let /x G C and co G s^ with a;(x) 7^ 0, < x < 1, be given. Put 



d u' d u + 1/2 u' 

dx Lo ax X oj 

d J 
d2:=--^ + —. 
ax oj 



(3.14) 



Then 



rii2 := aid2 = — -r-^ + 



dx"^ 



UJ 



rf2 /i2-l/4 ^/i + l/2w' u" 



, 2 + 2 + 2^^^^ + ^, (3.15) 

dX^ X^ X UJ LO 



H2i:=d2d, = -— + 2(^ 



d^ /a;'\2 cu" 



d' ^ (/. + !)- -1/4 ^ 2^i±V£E + 2fE)^ - ?, (3.16) 



dx^ x"^ X u \u J u 
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thus we have Hu = /^ + X"^Vi2, H21 = l^+i + X~^V2i and using uJ E s^ one directly 
checks V12, V21 G "^ (see Definition 1.2). 

Hence for a given V G fjicu the operator H = l^ + X^^V can be factorized in the 
form H = did2 with di, ^2 as above if and only if there is a solution to the homogeneous 
equation Hoj = with u{x) 7^ for < x < 1. Indeed, if u solves Hu = and 
is nowhere vanishing in (0,1], one directly verifies from the first line of (3.15) that 
H12 = did2 coincides with H. Conversely, given a factorization oi H = did2 with di,d2 
as in (3.14) it is immediate that d2UJ = and thus Hu = 0. 

By Proposition 2.6, oj is then of the form co{x) = x^'^"''^/^a3(x) with u G £/. 

The problem is that such a nowhere vanishing solution does not necessarily exist. 
However, we will be able to reduce the calculation of the zeta-determinant to the calcu- 
lation for factorizable operators. In fact, the main essence of the Proposition 3.5 below 
is that although H itself may not be factorizable, it becomes factorizable after adding 
a suitable L^^ (0, 1) potential. For such perturbations a variational formula for the 
zeta-determinant will be established subsequently. 

Next we investigate the separated boundary conditions for dj. For dj we can choose 
four possibly different closed extensions with separated boundary conditions: For p,q & 
{a,r} denote by dj^pq the closed extension of dj with boundary condition p at the 
left end point and boundary condition q at the right end point. Here, r stands for 
the relative boundary condition and a for the absolute boundary condition. More 
concretely, dj^rr = c^j.min is the closure of dj on C^(0, 1), dj^aa = c^j.max = {dj,min)* is the 
maximal extension. The domains of the mixed extensions can be characterized by 

nd,,ar) = {/ G ^(rf,,^ax) | /(I) = O}, 

^{dj,ra) = {/ G ^(rf,,^ax) | f (x) = 0{^ \ \og{x)\^l'') aS t ^ 0+}. 

For each choice of a closed extension D\ of d\ with boundary condition of the form 
aa^ rr, ar, ra we choose D2 to be the closed extension with dual boundary condition, i.e. 
rr, aa, ra, ar, for d2. If u is real then this means that D2 = D\. We summarize case by 
case the corresponding boundary conditions for if 12,-^21: 

Note that cu is a solution to the homogeneous differential equation if 125' = 0. In the 
notation of Section 1.4 we assume that 

uj = cos'i^o ■ 91 — sin-i^o ■ 92, (3.18) 

with < -i^o < TT, thus i?o,i?o'^ — 0- Moreover, we assume that Bi^^^^u = for a 
< -i?! < TT. We have to exclude -i?! = because in that case u{l) = and hence there 
would be a regular singularity also at the right end point. But see Remark 5.6. 
For future reference we now list the 

3.2. Separated boundary conditions for the factorized operator D1D2. 

Case I: Di = di^rr-,D2 = d2,aa- f G S!{DiD2) if and only if / G ^((i2,max) and ^2/ G 
■^(c^i.min)- Thus one checks that D1D2 = ifi2('*?0) ''^i) and i3'2i^i = if2i(0,0) is the 
Friedrichs extension of ii2i- 

Case II: Di = di^ra, D2 = '^2,ar- Then 

D,D2 = Hui'&O, 0), D2D^ = ii2l(0, TT - ^1). 
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Case III: Di = di^ar, D2 = '^2,ra- Then 

Case IV: Di = di^aa, D2 = '^2,rr- Then 

D,D2 = Hi2i0, 0), D2Di = if2i(vr - ^0, vr - ^1). 
We summarize the previous considerations in the following 

Proposition 3.5. Let V G i^^cu be given and let H = l^ + X^^V he the corresponding 
regular singular Sturm-Liouville operator. Suppose that we are given admissible sepa- 
rated boundary conditions Bjg.,j = 0,1, for H. Then for < -i^i < vr there exists a 
function uj{x) = x^^^l'^uj such that to E ^ (see Def. 1.2^, uj{x) 7^ for < a; < 1 and 
Bq 00^ = 0; Bi .Q^u = 0. Moreover, Hu = in a neighborhood ofO and in a neighborhood 
ofl. 

If 61 = we choose any < ?9i < vr and if < 61 < n we let 'di = 61. 

Putting Di = di^ra-, D2 = C?2,ar if ^ < ^q < It ,9i = 0; Di = di^rr, D2 = C?2,aa if 
< 60 < TT,6i > 0; Di = di^aa, D2 = d2,rr if Oq = 61 = and Di = di^ar, D2 = d2,ra if 

00 = 0,9i > we have 

H{eo,ei) = D^D2 + W (3.19) 

withWeLl^^{0,l). 

If V is real then u can be chosen to be real. 

Proof. We can certainly find an e > and solutions Uj,j = 0,1, to the differential 
equation Hg = on the intervals (0, e) resp. {l—e, 1] such that ojq^x) 7^ for < a; < e, 
a;i(x) 7^ for 1 — e < X < 1 and -Bo.So'^o = 0, Bi^^-^ui = 0. If V^ is real we may choose 
Uj to be positive. In any case we may choose a nowhere vanishing extension u to the 
whole interval with the claimed regularity properties. 

By construction D1D2 has the same boundary conditions as II{6o, 61) and there are 
neighborhoods of 0, 1 resp. on which the potential of D1D2 coincides with that of H, 
whence (3.19). D 



3.3. Comparison of Wronskians. For a factorizable operator H = D1D2 and given 
admissible boundary conditions we are now able to compare the Wronskians of normal- 
ized fundamental solutions of -D1-D2 and D2D1. 

Proposition 3.6. Let Rei^ > and let uj{x) = x~''^^^'^co{x) with to G £/ , co{x) ^ for 
< a; < 1. Put di,d2 as in (3.14) with fi = —v. Assume that Bq^^^^u = 0,Bi^^-^^u = 
with admissible boundary conditions Bj^^., < -^j < it for H12 = ^1^2- Choose Di,D2 
as in Proposition 3.5 with 9q = 'do and 61 =0 or 61 = "di. Let ifz,'ipz be normalized 
solutions for {D2D1 + z)g = 0. 

Case I: < 'di = 61 < TT . Then ip^ := ^^^diip, ip^ = —diipz are normalized solutions 
for {D1D2 + z)g = 0. Furthermore, we have 

W{^z, ^z) = ^^^(V^- ^^)- (3-20) 

D2D1 is invertible and the kernel of -D1-D2 is one-dimensional, hence spec D1D2 = 
specL'2^1 U{0}. 
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Case II: 9i = 0. Then (p^ := (2 — 2h')~^di(f, ipz = —z~^di'4'z ore normalized solutions 
for (-D1-D2 + z)g = 0. Furthermore, we have 

W{^,, ^,) = j^W{tP,, V,). (3.21) 

D2Di,DiD2 are both invertible and specD2Di = specDiD2. 

Proof. In view of (3.15), (3.16) the characteristic exponents of did2 are ±z/ + 1/2 and 
the characteristic exponents oi d2di are ±(— i/ + l) + l/2. diip^ satisfies {did2 + z)di(pz = 
di{d2d\ + z)^pz = 0. Furthermore, since (fz is normalized at for -D2-D1 + ,2 we have in 
the notation of (1.29) v?2(x) ~ x^^'^~'^ as x — )■ and using (3.14) we obtain 

diipzix) ~ (3/2 - z/ + 1/2 - u)x'/^-'' = (2 - 2u)x-''^'/\ (3.22) 

hence 2^^^^^^ is normalized at for -D1-D2 + -z as claimed. This applies to both cases 
I and II. 

Now consider diipz which also solves {did2 + z)diipz = 0. 

Case I: < -^i = 61 < IT. Then D2D1 = H2i{0, 0) and hence ipz being normalized at 
1 means ^^^(1) = 0,^^(1) = -1. Then ^1^^2(1) = ip'^i'^) + (^^^)(1) = "1 and hence 
—ditpz is normalized at 1. 

We now find for the Wronskian 

W{i)z,Vz) = ——--W{di^z,di(fz) = ———{di^z{diifz)' - {diiJz)'diipz) 

1 —z 

= {di'ilJz{.d2diipz) - {d2diiJz)di<^z) = {{.dii)z)'{>z -'^Pzdi'^z) 

= ^^W^(V^.,^.)- (3.23) 

Case II: 61 = 0. Then D2D1 = i72i(0,7r - ^i) with < n - ^i < tc. Thus ipz 
being normalized at 1 means ^2(1) = 1 and Bi^T^^^^ipz = 0. Then ipz G &{Di = di^ra) 
thus = diipzi^)- Hence (diipzYi^) = — (<^2C?i'0z)(l) = zipzi'^) = z. Thus —z^^diipz is 
normalized at 1 for H12. Now the same calculation as in (3.23) yields (3.21) and the 
Proposition is proved. D 



4. The asymptotic expansion of the resolvent trace 

Standing assumptions. Let Z^ := — ^ + '^ ~2 be the regular singular model oper- 
ator. In this Section we assume z/ to be real and non-negative. 

4.1. The Dirichlet condition at 0. Let Bi^e = sin^-/'(l) + cos^-/(l) be a boundary 
operator for the right end point and let Li, = L^{0, 0) be the closed extension of ly with 
domain 

^(L,) = {/ G ^(/,,^ax) I C2(/) = 0, 5i,,/ = 0}. (4.1) 

The following perturbation result will be crucial for establishing the resolvent trace 
expansion Theorem 1.4 (cf. [Les98, Lemma 3.1]). 
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Proposition 4.1. Let u > 0. Then L^ is self-adjoint and hounded below. 
Let W be a measurable function on [0, 1] such that 



(4.2) 



W^eL'^[0,l], ifz/>0, 

WHog{^) e L^[0,1], ifz/ = 0. 

Then for Zq > niaxspec(— L^,) there is a constant C{zq) such that for z ^'RjZ > Zq, we 
have for the Hilbert- Schmidt norms 

f (- + \l''^ \W{x)\^dx + i \}, ^W{x)?dx\, if z/ > 0, 

1 (i + /o^^ \W{x)\'\ log{xz)\dx + \ jlf^ l\W{x)\Hx), if ^ = 0, 

=: R{z). 

If (4.2) is replaced by W^ G %, i.e. 

'iy2iog(_) e^i[0,l], ifz/>0, 
WHoi{^eL^[QM ifz/ = 0, 

then 

lim R{z) = 0, (4.5) 

^-\R{z)\dz < oo. (4.6) 



(4.4) 



z—^oo 

fOO 



20 ^ 



Proof. The boundary conditions for L^, are separated and admissible. Therefore, L^, is 
self-adjoint. We will see below that the resolvent is Hilbert-Schmidt. Thus L^, has a 
pure point spectrum. An eigenfunction satisfying L^/ = A^/, A e MUiM, is therefore a 
multiple of y/xJ^{Xx), where J^, denotes the Bessel function of order u [Wat95]. From 
the known asymptotic behavior of the Bessel functions with imaginary argument one 
deduces that L^, has at most finitely many negative eigenvalues and hence is bounded 
below. 

The kernel ki,{x,y; z) of the resolvent {L^ + 2^)~^ is given in terms of the modified 
Bessel functions I^, K^ 

ku{x, y; z) = y/xyl^{xz) \K^{yz) - f3{z)I^{yz)j , x <y, (4.7) 

where f3{z) is determined by the requirement Bi^ek^x,-; z) = (cf. [BrSe85]). One 
finds 



m 



(cos 6^ + i sin 9)K^{z) + zKl{z) sin 9 
(cos 6* + I sin e)I^{z) + zl'^^z) sin 9 
(cos6' + {\ + v)sm9)Ky{z) - zK„+i{z)sm9 



(4.8) 



(cos6' + (| + z^) sin6')/;^(z) + zlyj^i{z) sin 6* 
where in the last equation we used the recursion relations [Wat95, 3.71] 

zll{z) = zh+i{z) + ul,{z), zKl{z) = -zK,+i + uK,{z). (4.9) 
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Recall the following asymptotic relations for the modified Bessel functions [Wat95, 
7.23] 

^pf z ^ oo, (4.10) 

and [Wat95, Sec. 3.7] 

1 „ , , (2''-^T(u)z-'', if z/^0, , , 

2*^1 (z/ + Ij I— logz, if z/ = 0, 

for the notation ~ see (1.29). From the asymptotics as z — )■ cxd one infers 

f3{z) = 0{e-^'), z-^oo. (4.12) 

To prove the estimate (4.3) we fix zq > maxspec(— L,^) and find for z > zq 

x-^\W{x)\'^K{x,x;z)dx (4-13) 



1 .1 

\W{x)\^I^{xz)K^{xz)dx - /3{z) / \W{x)\^\I^{xz)\^dx. 

Jo 

In the following C denotes a generic constant depending only on zq and u. We split the 
integrals into an integration from to 1/z and from 1/z to 1. In the first regime (4.11) 
yields 



\lJxz)K,{xz)\ < {^] , ,, '^''^O' (4.14) 

' ^ ^ ^ ^' - V|log(x2)|, if z/ = 0, ^ ^ 



and |/,y(x2;)| < C. Thus, 

-1/z 



X ^\W{x)\^ky{x^x\z)dx 



-1/z 



< 



C j^^"\W{x)\^dx, ifi/^0, 



(4.15) 



-1/z 



C j''"\W{x)\'^\\og{xz)\dx, ifi/ = 0. 



For 1/^ < X < 1 we apply (4.10) and find \Iy{xz)Ky{xz)\ < £, \Iu{,xz)\'^ < ^e^^^ 
Thus 



Ji/i 



,1 r' 1 



|iy(x)r|/^(x2)ir^(x2)|cia;<C- / -|W^(x)|'cia;, (4.16) 

/z Z Jl/z X 
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and in view of (4.12) 

1/3(^)1 r \W{x)\'\I,{xz)\'dx 

Jl/z 

<Ce-^' / \Wix)\^—e^'''dx 

<c(e-'- [ -\W{x)\^dx+- [ \W{x)\^dx' 

\ Z Jii-, X Z Ji/2 

< C-, Z > Zq. 

z 

The estimate (4.3) now follows from (4.13), (4.15), (4.16), (4.17). 

Under the assumptions (4.4) we apply Lemma 2.4 to --R(-) since J°° -\R{z)\dz = 

lo^'' ^|i?(^)Mx and conclude (4.5), (4.6). D 

We return to the discussion of the operator H = l^ + X~^V; recall that X denotes 
the function X{x) = x. We have seen in the previous Proposition that Li, is a bounded 
below self-adjoint operator. In fact L^ is the Friedrichs extension of ly restricted to the 
domain 

^{Q = {/ e ^0^(0, 1] I B,,ef = 0}. (4.18) 

We now want to construct the Friedrichs extension of H on ^{l^) and compare its 
resolvent to that of L^; cf. [Kat95, VI, 2.3]. The problem is that the domains of Li, 
and of the Friedrichs extension of H on !^{lu) are not necessarily equal. This is because 
the domain of L^ contains functions f{x) with f{x) ~ x'^^^^'^ as x — )■ 0. For such a 
function, X~^Vf is not necessarily in L^[0, 1]. 

Proposition 4.2. Let H = l^ + X'^V with V e% (cf. Def. 1.2^. Moreover, let &{Q 
be given by (4.18), < ^ < tt, and let q{f,g) := {luf,g) be the form of the operator ly. 
Then the form 

v{f,g) := (X-V/,^)i2[o,i], f,g e m) (4.19) 

is q-bounded with arbitrarily small q-bound b. 

Proof. We compute for any g G ^{lu) and z > zq 

Hg,g)\=\\x-'/'\V\'/'g\\l. 

<\\x-'/'\V\'/'iLy + zY'/Y{iLy + z')g,g). 

Now Proposition 4.1 implies, that for any 6 < 1 there exists z G M_|_ sufficiently large, 
such that 

Hg,g)\ <b{{Ly + z^)g,g) = bz^\\g\\l2+bq{g,g). D 

The quadratic form q is bounded below and closable with closure Q. By the second 
representation theorem [Kat95, IV, 2.6 Theorem 2.23], we have 

^{Q) = ^{{Ly + ziy/'). (4.20) 
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As a consequence of Proposition 4.2 we find in view of [Kat95, VI, 1.6, Theorem 
1.33] that (g + f ) is a sectorial form with 

^(q + ^) = ^(Q) = ^{{L, + zlY^^). (4.21) 

By the first representation theorem ([Kat95, VI, 2.1, Theorem 2.1]) it determines 
uniquely a closed m-sectorial extension H{0, 6) oi H = l,^ + X~'^V, with domain given 
by 

= {fe nHra.^) I C2(/) = 0,i?i,,/ = 0}. 

Note that the functional C2 (as well as ci) depends on the potential and the C2 in (4.22) 
is the one associated to H. 

Theorem 4.3. The operator H{0, 9) is m-sectorial, in particular spec H{Q, 9) is a subset 
of a sector {^ G C | | arg(^ — 'n)\< a}, for some fixed angle a G (0, 7r/2) and 77 G M. Its 
resolvent is trace class and 



Ri{z) : = 


\iH{0,9) 


+ 


z) ^ -{L^ + z) ^\tr, zeR+,z>max{-r],0) 


(4.23) 


satisfies 






limzRi{z) = 0, zeR+, 

2— s>oo 

/ \Ri{z) dx < 00. 

Jzo 


(4.24) 
(4.25) 


Furthermore 


J 









-1 m / r , \-l 



Ti{H{0,9) + z) ' = Ti{L, + z) +R2{z 



a b ^, , (4.26) 

'z z 



where 

^ = 1 b = -hu + ^^,{B^,e)) = H[" ^ \\ '!^'J°' (cf. (1.28)) (4.27) 

2 2^ l~2(^~2)' II < t^l < TT. 

The remainders R2{z),R^{z) satisfy (4.24) and (4.25) and therefore the zeta- 
determinant of H{0,9) is well-defined by the formula (see (1.12) and Figure 1). 

\ogdet^H{0,9) = -JTT{{H{0,9) + z)-^)dz. (4.28) 

Proof. The operator H{0, 9) is m-sectorial, as it arises from the sectorial form {q + v), 
see [Kat95, VI. 2, Theorem 2.1]. Since by Proposition 4.1 we have 

\\m\\x-'/'\V\'/'iL^ + z)-'/'\\us = 0, 
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we may invoke the Neumann series to obtain 

iHio,e) + z)-'-iL, + zY' 

= 5^(-l)"(L, + z)-"^ [{L, + z)-^^x-'V{L, + zy'^YiL, + z)-"^. (4.29) 

n>l 

There is a httle subtlety here since &(^H{0, 9)) does not necessarily equal ^{Ly). How- 
ever, by Proposition 4.2 the forms of -ff (0, 9) and Ly have the same domain. This is 
used decisively by writing {L,, + 2;^)"^/^ at the beginning and at the end of (4.29). 
We estimate the trace norm of the individual summands by 

||(L, + ^)-i/2f . ||(L, + z)-^'^x-'V{L, + z)-^'m 

< ||(L. + z)-^\\x-'/'\V\'/\U + .)-^/^||Ss ■ UU + zr^'x-'^'lVnins 

< Cz-^R{zY, 

where R{z) = \\x'^l^\V\^'''{L, + z)-^/'||hs ■ II (^. + z)-^'^x-^'^\V\^'''\\ns- The claim 
about Ri{z) now follows from Proposition 4.1. 

The first line of (4.26) follows since |i?2(2;)| < Ri{z). As for the second line of (4.26) 
we note that Tr(Li, + z)~^ has a complete asymptotic expansion as 2; — )■ 00 [BrSe85], 
in particular 

Ti{K + z)-' = ^ + - + 0{z~^/Hogz), 

^ ' 'z z 



with a, h as in (4.27). 

For the claim about the zeta-determinant see Section 1.2. D 

4.2. General boundary conditions. We now extend Theorem 4.3 to general bound- 
ary conditions at 0. Recall that is in the limit point case if and only if z/ > 1. 
So the following discussion is of relevance only in the case v < \. The case z/ = 
bears more difficulties (see [FMPS03], [KLP06]) and therefore we assume from now 
on < z/(< 1). The difficulty then is that for < 6'o < vr the resolvent of lu{9o,9i) 
does not absorb negative x powers as the operator 1^(0, 9i) does. Therefore, we do not 
have (4.3) at our disposal and hence the resolvent of H{9o, 9i) cannot be constructed 
as a perturbation of the resolvent of lu{9o, 9i). Instead we will employ the results about 
factorizable operators in Section 3.1. However we have to impose a slight restriction on 
the class of potentials: 

Definition 4.4. Let V & 'Vy and let H = l^ + X~^V be the corresponding regular 
singular Sturm-Liouville operator. V is called of determinant class if for any pair of 
admissible boundary conditions Bjg the operator H{9q, 9i) satisfies for z > zq^z & R+, 

\\{H{9,,9^) + zY^\\=0{\z\-^), (4.30) 

\\{H{9,,9,) + z)-\ = 0{\z\-^"), (4.31) 
and for any (/) G C^(0, 1] 

\\v{H{9o,9^) + z)-'\\^,^^,=0{\z\-''^). (4.32) 
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Here, || ■ ||i2_>j;^i denotes the norm of a map from L^[0, 1] into the first Sobolev space 
H^[0, 1]. We denote the set of determinant class potentials by Y^'^^. 

We note some consequences and give some criteria for V being of determinant class. 

Lemma 4.5. Let V e f^^"^ and let W e Lcomp(0, 1] with supp W c[6,l],6>0. Then 

\\WiH{eo,ei) + z)-'\\ <Cs\\W\\l2 \z\-^/\ z>Zo. (4.33) 

ForW e L°°[0,1] we have 

\\W{Hi9o,9i) + z)-'\\<C\\W\\oo\z\-\ z>zo. (4.34) 

Proof. Choose a cut-off function (f G C^(0, 1] with <f{x) = 1 for x > 6. Then (4.30), 

(4.32) and the complex interpolation method [Tay96, Sec. 4.2] yield for < s < 1 

\\^{Hi9o,9,) + z)-'\\^,^^^<Cs\z\-'^'/'. (4.35) 

By Sobolev embedding we have -^''[0, 1] C C[0, 1] for s > 1/2 and thus for these s mul- 
tiplication by W is continuous H^ — )■ L^ with norm bounded by Cs||iy||L2. Combining 
this with (4.35) gives 

\\W{H{9o,9,) + zy'W^,^^, < CsAWWl^-'^'^^. (4.36) 

(4.33) follows by putting s = 2/3, (4.34) is obvious from (4.30). D 
Lemma 4.6. Let V e yj^""'. If W = Wi + W2, Wi e L~[0, 1], 1^2 e Lcomp(0, 1] then 

v + xw e r^^y ^ 

Consequently, ifV\ G 'y^'^^ ■, V2 & "Vu and Vi{x) = V2{x) for almost all x in a neighbor- 
hood 0/0 then V2 £ 'f'^'^^ ■ Furthermore, there is a constant depending only on H{9q,9i) 
and the support of W2 such that for z > zq 

\\{H{9,,9,) + W + z)-' - {H{9o,9,) + z)-\ 

<c{\\w,\\^ + \\w2\\l^) i^r'/^ 

Proof. It follows from Lemma 4.5 that for z large enough we can employ the Neumann 
series 

{H{9oA) + W + z)-' - {H{9o, 9{) + z)-^ 

= Y,{-lT{H{9,,9,) + z)-\w{H{9,,9,) + z)-'Y ^^'^^^ 

n=l 

and (4.30), (4.31), (4.32) follow for H{9o,9i) + W; also (4.37) immediately follows. 
The second claim follows from the first with W = X~^(V2 - Vi) e L^^^^^O, 1]. D 



Proposition 4.7. Let V G '% be real valued in a neighborhood ofO. Then V E 'fl 



dct 



Together with Lemma 4.6 this shows that at least potentials of the form V + X, where 
V E y^ ^s real valued and A G C, are of determinant class. 



^Note that then H + W = l^ + X~\V + XW). 
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Proof. In view of Lemma 4.6 and Proposition 3.5 we may change V outside a neighbor- 
hood of such that V becomes real valued everywhere and such that H{9o,9i) = D*D, 
where D is a closed extension oi d = —-^ + uj'/uj. For the properties of a; see Proposition 
3.5. Note that since V is real valued we may choose uo to be real valued, too and hence, 
in the notation of Proposition 3.5, Di = D^- 

Since D*D is self-adjoint, elliptic and non-negative (4.30), (4.32) follow immediately 
from the Spectral Theorem. If ^o = then (4.31) follows from Theorem 4.3. If ^o 7^ 
then by Proposition 3.5 the operator DD* has Dirichlet boundary condition at 0. Hence 
by Theorem 4.3 the estimate (4.31) holds for DD* . Since for a non-negative operator the 
estimate (4.31) depends only on the spectrum and since specDD* U {0} = specD*D U 
{0} we reach the conclusion. D 

Next we prove two comparison results for the asymptotics of the resolvent in the trace 
norm. These will then lead to an asymptotic expansion of the trace of the resolvent 
for H{9o,0i) for arbitrary admissible boundary conditions and all determinant class 
potentials. The technique used in the first comparison result is well-known for elliptic 
operators with smooth coefficients on manifolds (cf. e.g. [LMP09, Appendix B]). We 
have to be slightly more careful here due to the low regularity assumptions on the 
potential. 



Proposition 4.8. Let Vj G yj^''\j = 1,2 with V2 ~ Vi G 11^^^(0,1], that is there 
is a 6 > such that Vi{x) = V2{x) for < x < 6. Let Hj = ly + X~^Vj he the 
corresponding regular singular Sturm- Liouville operators and let Bg^^Bg^ resp. Bg , he 
admissihle houndary conditions for Hj. Then there is a zq > such for any 6' < 6 and 
z > zq the difference [Hi{9o,9i) + z) — (^H2{9o,9i) + z^ restricted to L'^[0,6'] is of 
trace class and 



H,{9o, 9i) + z)'' - {H2i9o, 9i) + z)-') 



L2[0,5']lltr 



0{\z\ 



-3/2N 



z > zo,z e 



Proof. We choose cut-off functions (^, ip 
identically one over [0,6'] and 

• supp(0) C supp (■?/'), 

• supp((/)) D snpp (dip) = 0. 



C^[0,S), cf. Figure 2, such that they are 




Figure 2. The cutoff functions (p and ip. 
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In particular these conditions yield tpcj) = (f). In this proof we will write for brevity 
Hi instead of Hi{9o, 9i) and H2 instead of H2{0o, 9i). 
We now consider 

R{z) := ^ [{Hi + z)-' - {H2 + z)~']<j). 

R{z) maps into the domain of Hj and on the support of ip the differential expressions 
Hi and H2 coincide; moreover ip^iHi) = ^'^(ifa)- Thus {Hi + z)R{z) = [iJi,^]((iJi + 
z)~^ — (H2 + z)~^y Arguing similarly for R{z)* and taking adjoints one then finds 

{Hi + z)R{z){H2 + z) = [-dli,]{{Hi + z)-' - {H2 + z)-')[dlAi 

where [■, ■] denotes the commutator of the corresponding operators and any function is 
viewed as a multiplication operator. Hence 

R{z) = {Hi + z)-'[-dlM{{Hi + zr'-{H2 + z)-')[dlct>]{H2 + z)-' 

and thus 

l|i?(z)llt.<||(i/i + ^)-^llt.(||[9^,^](i/i + z)-^ll + l|[a^,^](if2 + ^)-^ll)- 

■\\[dl<P]{H2 + z)-\ 

By (4.31) we have \\{Hi + z)-^\\i, = 0{\z\-^/^). Let / denote ip or 0. Then [d^ f] is 
a first order differential operator whose coefficients are compactly supported in (0, 1), 
hence it maps H^[0, 1] continuously into -^comp(0' !)• Therefore by (4.32), with a cut-off 
function x ^ C^{0, 1) with % = 1 in a neighborhood of supp([9^, /]), 

\\[dlj]{H, + zr'\\ < \\[dlJ]\\H^^L4xiH, + zy'u.^j,. = oi\z\-'/'^ 

ioT z>zo. Hence \\[dl,(f)]{Hj + z)-^\\ = 0{\z\-^/^) and \\[dl,tlj]{Hj + z)-^\\ = 0{\z\-^/^) 
and the proposition is proved. D 

We note that in this proof the estimate (4.31) was used only for Hi. 

Proposition 4.9. Let Vj G r;^°S Hj = l^ + X'^Vj,j = 1,2, and let B^^^, Bq^Oo, Bi^e^ 
be admissible boundary conditions. Then for any 6 > 

II {{Hi{9o, 9i) + zr' - {H2(9o, 9i) + zY^) \^,y,^,^\l^ = 0{\z\-^'^), z>zo,ze M+. 
Proof. Fix (5 > and put 

c/2 /Z/2-1/4 1 



H. 



+ a^-^^ + ^^1) l[V2,i] =: A + 9, (4.39) 



withgGL2„^p(0,l],A:=-^. 

Exactly as in Proposition 4.8 one now shows 

II {{Hi{9,, 9i) + z)-' - {H,{9„ 9i) + z)-') |^,j,_ J^^ = 0{\z\-'"), (4.40) 

for z > Zq,z & ]R+. Furthermore, an elementary calculation involving the explicitly 

computable resolvent kernel of A(6'o,6'i) shows ||g(A(^05^i) + -2)~^||tr = 0{\z\''^/'^). A 
Neumann series argument then gives 

II {{H,{9o, 9i) + zY' - (A(^o, ^1) + ^V)\l2^saL = Oi\zr'^')- (4.41) 
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(4.40) and (4.41) imply ior z > zo,z eR+, 

II {{H,{eo, e,) + z)-' - (A(^o, ^i) + ^)"') L.[5,i] lltr = 0{W'). (4.42) 

The same line of reasoning applies to H2 and hence (4.42) also holds with H2{6o,6i) 
instead of Hi{6o, 61), whence the result. D 

Theorem 4.10. Let V E Yif'^^,u > 0, and let H = l^ + X~^V be the corresponding 
regular singular Sturm- Liouville operator. Let < 6j < 71 (Oq = if u > 1). 

Then the resolvent of H{6q,6i) is trace class. Moreover, there is a zq > such that 
H{9o,9i) + z is invertible for z > zq and 

TT{H{eo, ei) + zy' = -^ + - + Rs{z), z > zo, 

where a 



z z 



2' 



b = -7;{f^{Bo,eo) + KBi,0,)) = < 



f-\u-\, if ^0 = ^1 = 0, 

-iz/ + i, if ^0 = 0,0 <^i <7r, 

\y-\, if 0<eo <vr,ei = 0, 

^\v+\, if 0< 00,^1 <vr, 



(4.43) 



(cf. (1.28)j is independent ofV, and the remainder Rs{z) satisfies 

lim zRsiz) = 0, ze M+, (4.44) 

/•oo 

/ \Rs{z)\dx < 00, zo> raaxspec{-H{9o,9i)) nR. (4.45) 

J Z() 

In particular the zeta- determinant of H {60, 61) is well-defined by the formula (1.12), 

logdet^L = - -/ Tr(/7(0, 9) + z)-^)dz. (4.46) 

Proof. By Proposition 3.5 we may choose a factorizable operator Hi{9o,9i) = -D1-D2 
such that there is a 5 > such that the coefficients of Hi and H coincide on the 
interval [0,6]. Here, Di,D2 are appropriate closed extensions of the operators di,d2 in 
(3.14) with fi = —V. Then by Propositions 4.8, 4.9 we find 



(i7(^o, ^1) + z)-^ - {D1D2 + z)-^) IL^ = 0{\z\-'"^), z>zo,zE M+, (4.47) 

and hence 

Tr(i7(0o,^i) + zY^ = Tt{DiD2 + z)'' + 0{\z\-''/^). (4.48) 

We now have to discuss the four possible cases listed in Section 3.2, see also (3.15), 
(3.16): 
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Case I: Di = di^rr, -D2 = d2,aa- Then -D1-D2 has a one-dimensional null space and -D2-D1 
is invertible. Applying Theorem 4.3 to -D2-D1 we obtain 

Tr{DiD2 + z)~' = Ti{D2Di + z)'' + z'^ = ^ + ^ + R^{z), (4.49) 



where Rz{z) has the claimed properties (4.44) and (4.45) and a = 1/2, b = —1/2(1 — 
1^ + 1/2 — 2) = l/2(z/+l/2). Note that in formulas involving -D2-D1, according to (3.16), 
the u has to be replaced by 1 — z/. 

Case II: Di = di^ra,D2 = '^2,ar- Then D1D2 and D2D1 are both invertible and hence 
Tt{DiD2 + z)^^ = Ti{D2Di + z)^^, and we can proceed as in Case I. 

In the remaining cases III {Di = di^ar, D2 = '^2,ra) and IV {Di = di^aa, D2 = '^2,rr) 
one can apply Theorem 4.3 directly to D1D2. □ 

5. Variation of the regular singular potential 

In this section we discuss the behavior of the fundamental system of solutions un- 
der a certain variation of the potential and derive a variational formula for the zeta- 
determinant. 

Standing assumptions. Let z/ > 0, 1/ G f^^^^ and let Wr, G L°°[0, 1]+Ll^^p{0, 1] be a 
family of functions depending on a real or complex parameter rj. To avoid unnecessary 
technicalities we assume that Wr, is of the form Wr, = Wi r, + W2 n where Wi ^ G 
L-[0,1], 1^2,, eL2„^p(0,l] satisfy 

(i) rj I—)- Wi^r) is differentiable as a map into the Banach space L°°[0, 1], 
(ii) there is a fixed 5 > such that supp W2,n C [5, 1] and i] 1— )■ W2,r) \ ,^ ^, is differen- 
tiable as a map into the Banach space -^^[5, 1]. 
For notational convenience we assume Wq = and put V^i := V + XW^^ and 

d^ 
Hr, := /. + X-% = l, + X-\V + XW^) =: Hq + W^ =: -— + g^. (5.1) 

rjo = serves as a base point for a perturbative construction of a fundamental system. 

5.1. Fundamental solutions and their asymptotics. According to Theorem 2.1 

let gi^r] be the unique solution of the ODE Hr^gir, = with gi,n{x) ~ x^^, as a; — )■ 0+. 
Note that the second solution (yf2,r? in Theorem 2.1 is not uniquely determined by the 
requirement g2,ri{x) = —^x'^^, cf. Remark 3.2. Since the solutions now depend on the 
parameter rj, the choice of g2,n becomes important. Before we specify g2,Tj we discuss 
the dependence of gi^r^ on rj. To do so recall the operator Ki, from (2.11) in Section 2. 
For a > consider the Banach space X°'C[0, 1] with norm 

||/|U:= sup |a;-"/(x)|, (5.2) 

0<i<l 

and the Banach space C^[0, 1] consisting of those functions in / G (X°C[0, 1]) nC^(0, 1] 
with /' G X"-iC[0, 1] and norm 

ll/llc^:=ll/IU+ll/'IU-i- (5.3) 
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For / e X°C[0, 1] the inequality (2.15) gives 

\iKVrfix)\ < x'^T^WfUr \V{y)\dyy, v ^ 0, (5.4) 

V'rnl \Jq ) 

thus KyV is a bounded operator on X"C[0, 1] with spectral radius zero. Furthermore, 
for / e X"C[0, 1] (cf. (2.19)) 

\[KMf){x)\<x'^-\\no.f \V{v)\dy. (5.5) 

hence K^V maps X"C[0, 1] continuously into C^[0, 1]. 
Recall from (2.17) that g\^n{x) = x''"^{l + 0,,(a;)) with 

0, = (/ - KV,)-'KV,1. (5.6) 

Consequently 0^ is different iable in rj and 

d,ct>, = {I-K,V,)-'K{Xd,W,)l 

+ (/ - KV,)-'K,{Xd,W,){I - KV,)-'K,V,1. 

Since d^^Wrj is bounded near 0, the operator K^{XdrjWrj) maps C^[0, 1] continuously 
into C^+2[0, 1] and hence we have proved 

Lemma 5.1. Under the assumptions stated at the beginning of this section, gi^^i is dif- 
ferentiable in rj with drjgi^rj{x) = 0{x'^^'^'^),dr]g'iJx) = 0{x'^^^^) as x ^^ 0+. Moreover, 
the 0-constants are locally uniform in t] and hence gi^ni^) — gi^rioi^) = 0{x'^^'^'^). 

After these preparations we can discuss the second fundamental solution g2^n- For rj 
in a neighborhood of we can fix xq G (0, 1) such that gi^rji^) 7^ for < x < Xq. For 
these X we note 

giA^y^ - 9i,oixy^ 

_ [gi,^{x) + gi^o{x)][gi^o{x) - gi^^jx)] _ ^_^ (5.8) 

{gr,,{x)ng,,,{x)f ^^"^ ^' '^^'+' 

where the 0-constant is independent of rj. Hence g^^ — g^^ is integrable over (0,a;o] 
and we put for x G (0, Xq) 

rxo 

g2fi{x) = gi,o{x) / gifiiyy^dy, (5.9) 

J X 

,{X) r[g,^^(y)-^-g,M'']dy+^-^g2,o{x). (5.10) 



Jo gi,o[x) 

From (5.8), (5.9) and (5.10) we immediately get 

Lemma 5.2. gi^n,g2.r] is a fundamental system of solutions for the ODE Hr^g = 

satisfying (2.5), (2.6). Moreover, g2^n is also differentiate in rj and we have for z/ > 

^2,,(x)=^2,o(x) + 0(x'^^+2)^ ^5_1^) 

d,g2,,{x) = 0{x^-^''), 5,^2,,(^) = Oix^-^^), (5.12) 

as X — )■ 0. For V = the estimates are 0(x^/^~''logx) = 0(x^/^logx), 0(x^2"'"^logx) = 
0(x^/^logx), ©(x'^^^^ logx) = 0(x^/^logx), respectively. 
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Lemma 5.1 and 5.2 imply: 

Corollary 5.3. For u > we have the following asymptotics for the Wronskians 
W{gj,^, d^gk,r,) = 9j,v9'k,^ ~ 9j,rj9k,v,j, k = 1,2, as x ^ 0+: 

W{g,,„d,g,,,){x) = 0{x'''^+'), (5.13) 

W{g2,,,d,g,,,){x) = 0{x'^+''^+') = 0{x'), (5.14) 

W{g,,,,d,g2,,){x) = 0(x^^+^^+i) = 0{x'), (5.15) 

W{g2,,,d,g2,,){x) = 0{x^'''+'). (5.16) 

If u = then the estimates are 0{x^logx) in all four cases. 
Hence for u > and all j,k = 1,2, we have 

\imW{gj^rj,dr,gk,^)ix) = 0. 

Now we are ready to state the variational result which generalizes [Les98, Prop. 3.4] 
to arbitrary boundary conditions and to more general potentials: 

Theorem 5.4. 1. LetO < u <1, V e YJ^''^ and let r] ^ Wr, e L~[0, 1] + i^comp(0' 1] 
be differentiahle in the sense described at the beginning of this section. Furthermore, 
let < 9j < TT,j = 0, 1 and let H^ = l^ + X~^V + Wr,. Fix tjq and let gj^,^ be the 
fundamental system constructed above, relative to the base point tjq (gj,r,o plays the role 
of the gjfi above). 

Then we have H^{6q,6i) = Hrn^(9Q,9i) + Wr, — Wr,(^. Moreover, if Hn^{9Q,9i) is in- 
vertible then t] ^^ logdet^ if^(0o7 ^i) is differentiable at tjq and if (pr^^ipr, denotes a fun- 
damental system which is normalized for the boundary conditions Bj^q.,j = 0, 1 we 
have 

^\^\ogdet^H,{9,,9,) = ^IjogWiij,,^,). (5.17) 

2. Let u > and let r] ^-^ Vri & i^u be differentiable (recall from Def. 1.2 that 1^ is 
naturally a Frechet space). Let Hj^ = l^ + X'^Vj^, < 9 < n. If -ffr;o(0, ^o) is invertible 
then rj i— )■ logdet^-ff^(0, ^) is differentiable at tjq and formula (5.17) holds accordingly. 

Proof. 1. Let < z/ < 1. By Lemma 5.1 and Lemma 5.2 we have 



9iJx) -gi„,,{x) = 0(x^/^ 



a;^0. (5.18) 

92,7j{x) - g2,no{x) = 0{x ' ), 

Hence by Theorem 3.1 the domain of -ffr,,max as well as the functionals Ci,C2 are inde- 
pendent of rj. Thus we have indeed Hr^{9o,9i) = i?^^(,(6'o, 6*1) + Wff — Wr^y The proof 
of Lemma 4.5 shows that Wr, is i7^o(6'o, 6'i)-bounded and the assumptions on the map 
r] H^ Wn then imply that rj i->- H^{9q, 9i) is a graph continuous family of self-adjoint op- 
erators; in particular there is an e > such that Hn{9Q, 9i) is invertible for \f] — f]Q\ < s. 
From now on we assume |?7 — r^ol < £. 

From the estimate (4.37) we conclude that 

log det^Hr,{9Q, 9i) - log det^ Hr,, {9o, 9{) 

= -fTT{{H,{9o,9,) + z)-' - {H,,{9o,9,) + z)-')dz ^^'^^^ 
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where the integrand on the right is absolutely summable as it is 0{\z\^'^^^), z — )■ oo. 
Furthermore, according to our assumptions on VF^ we have 



d 

= -iHr,{9o,9i) + zy\^r,Wr,){H,,{eo,9^) + z)-\ 



■ [{H^iOo, di) + z)-^ - {H,, (00, 0i) + z'-' 



By (4.37) the trace norm of the right hand side is 0(|z|^^/^) where the O— constant 
is locally independent of t]. By the Dominated Convergence Theorem we may thus 
differentiate under the integral and find 

' TT{{H^{eo,ei) + z)-\dr,w^){H^{eo,ei) + z)'')dz 

r (5.21) 

^TT{{^,W,)iH,i9o,9^) + z)-')dz 
p cLz 

= TT{{dr,W^)H^{9o,9,)-'). 

Having established this identity we can now proceed as in the proof of [Les98, Prop. 
3.4], making essential use of Corollary 5.3. 
The kernel Gn{x, y) of H^{9q, 9i) is given by 

Grjix, y) = W{i)^, ip^y^ipr,{x)ip^{y), x <y, (5.22) 

W{iprj, V^r;) 7^ siuce i?r;(^0) ^i) IS iuvertiblc by assumption. Since gj^rj are differentiable 
in 7] (Lemma 5.1 and Lemma 5.2) so are ipni'4'v -'■^ fact, the normalization condition 
implies 

_ |-cot6'o ■5'i,r, + fi'2,,7, ifO<6'o<7r, 
'^'"1^1,,, ifeo = 0. (5.23) 

i^n = (^v9i,v + bv92,v, 

where a^, 6^ depend differentiably on rj. Differentiating the formula (pg^^ = q-q^e,r) with 
respect to rj gives 

dr^ip'^ = {dnqr,)(Pr, + q-qdnifr^ = {dr,Wr,)iPr, + Qr^drjifr^, (5.24) 

and hence 

d d (5.25) 
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Thus we find 

_ f^ d 

= W{^r„^rj)~HW{tpr„d^(p^){l)- lilll W {ip^, drj(p^){x)) . 



By Corollary 5.3 we have 
On the other hand 



lim W(tp„,dr,ipr,){x) = 0. (5.26) 

X-5-0+ 



Wid,iJr,,^,)il) = 0, (5.27) 

since ■?/'^ is normalized with ipr]{^) = and ip'Jl) = —1 in case of Dirichlet boundary 
conditions and with V'r;(l) = 1 in case of generalized Neumann boundary conditions. 

Note that in contrast to [Les98], the proof of relation (5.26) requires a careful as- 
ymptotic analysis of the fundamental solutions as summarized in Corollary 5.3. 

In view of (5.26) and (5.27) we arrive at 



drj 



d d 



and the proof of 1. is complete. 

2. For the proof of 2. we only have to note that by Proposition 4.1 we can estimate 
the trace norm of {H^{0,9) + z)-\d^Vr,){Hr,{0,9) + z)-^ by C\z\-^R{z) where R{z) 
satisfies (4.6) and the constant C is locally independent of rj. Thus we conclude the 
variation formula (5.21). The remaining arguments are then completely analogous to 
the proof of 1. D 

Remark 5.5. One can also prove a variation formula for the dependence of the zeta- 
determinant on the boundary conditions 6o,6i- For the variation of 6i at the regular 
end this is standard, see e.g. [Les98, Prop. 3.6]. For the variation of 9q the proof is 
much more delicate. Due to our approach via factorizable operators the result is not 
needed and therefore omitted. However, the factorization method does not extend to 
matrix valued potentials in a straightforward way. So, if one would like to generalize 
the results of this paper to matrix valued potentials with regular singularities then one 
would probably need to establish a formula for the variation of the zeta-determinant 
under the variation of the boundary conditions at the singular end. 

5.2. Proof of the Main Theorem 1.5. We are now finally ready to prove the Main 
Theorem 1.5. As in [Les98, Sec. 4] we first note that (1.30) is obviously true ii H{6o, 6i) 
is not invertible. Furthermore, if </?(-, ;z), ?/'(■, 2;) denote the normalized solutions for 
H(9q,6i) + z it follows from Theorem 5.4 (surely, for V G Y^*^^ the family z \^V + zX 
satisfies the standing assumptions of the beginning of this section) that det^(iJ(6'o, 6*1) + 
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z) and W{ilj{-,z),ip{-,z)) are holomorphic functions in C with the same logarithmic 
derivative. Hence it suffices to prove the formula for H[6o, 6i) + z for one z E C 

Let us now first assume that ^o = 0, i.e. at the left end point we have the Dirichlet 
boundary condition. Except for the low regularity assumptions on the potential this 
case was treated in [Les98]. From loc. cit. we will only use the result that the formula 
(1.30) holds for 6*0 = and V{x) = xz, i.e. for the operator /,^(0,6'i) + z. To reduce 
the claim to this case we consider Vr, := rjV. By Proposition 4.1, L^ := /j,(0, 9i) is self- 
adjoint and bounded below and from (4.3) we infer that -ffr?(0, Oi) + z := L^ + riX~^V + z 
is invertible for < r/ < 1 and z > zq. Hence we may apply the variation result Theorem 
5.4, 2. and we are reduced to the case V = and thus to [Les98]. 

Next we consider the case < ^o < 1- As noted before this necessarily means u < 1, 
since for u > 1 the left end point is in the limit point case. The case z/ = is beyond 
the scope of this paper and so we assume < z/ < 1. By Proposition 3.5 we have 
H{9o, di) = D1D2 + W with W e Ll^^piO, 1] and Di = d^^ra, D2 = ^.ar if Oi = 0, and 
-Di = di^rr,D2 = c?2,aa if ^'i > 0. Putting Wrj = riW,0 < r] < 1, we infer from Lemma 
4.5 that D1D2 + r]W + z is invertible for < r^ < 1 and z > zq, thus invoking again 
the variation result Theorem 5.4, 1. we are reduced to prove the formula (1.30) for 
the operator D1D2 + z. Note that -D2-D1 has the Dirichlet boundary condition at 
and hence (1.30) holds for D2D1 + 2; by the ffist part of this proof. We now look at 
the cases already discussed in the proof of Proposition 3.6. We use the notation from 
loc. cit., in particular ip^, 4'z denote a pair of normalized solutions for {D2D1 + z)g = 
and (fz = j^diipz, ipz = —diipz the corresponding pair of normalized solutions for 
(D1D2 + z)g = 0. 

Denote by /io,/^i the invariants defined in (1.27), (1.28) of the boundary conditions 
for -D1-D2. Denote by /i' the corresponding invariants for -D2-D1. 

Case I: < ^1 < TT. The kernel of -D1-D2 is one-dimensional and -D2-D1 is invertible. 
We have /^o = —^,fi'i = —1/2. /Xq = 1 — i',fi'i = 1/2. Thus, using the proven formula 
(1.30) for -D2-D1 and Proposition 3.6 



det^{DiD2 + z) = zdet^{D2Di + z) 

ZTl 



2/^o+Mir(/i^ + i)r(;u; + 1) 



W{tpz,^z 



2(l-z.)7r ~ ^ (5.28) 

W{^IJz,Vz) 



2M+m+2r(/io + 2)r(/ii + 2) 
2Mo+mr(/io + i)r(/ii + i)^^'^^' '^ 
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Case II: 9i = 0. Here -D1-D2 and D2D1 are both invertible and we have /xq 
1/2. /io = 1 - i^,fi[ = -1/2. Thus 

det^{DiD2 + z) = det({D2Di + z) 



-Z/,/ii 



TT 



2'^o+M'ir(/i^ + i)r(^; + i 

2(l-z/)7r 



-W^(^2,<^2 



2Mo+Mir(yLio + 2)r(/ii 



Vr(V^2,</'2 



(5.29) 



TC 



2Mo+Mir(/io + i)r(/ii + i 



■w^(^.,^.). 



The proof is complete. 



D 



Remark 5.6. We conclude by mentioning that Theorem 1.5 can be extended to po- 
tentials with regular singularities at both end points (and otherwise having the same 
regularity properties as the class Y^^^). The formula (1.30) remains the same. For the 
proof one first employs the factorization method we used here to arrange that, say at 
the left end point, one has Dirichlet boundary conditions. For this boundary condition 
a variation formula for the variation of the singular potential was proved in [Les98, 
Prop. 3.7]. This variation formula is still valid for our class of potentials and it allows 
to deform the parameter i/ to z/ = 1/2. Now one is basically in the situation with one 
regular end point and one singular end point and Theorem 1.5 can be applied. The 
details are left to the reader. 
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